Formation of u Centauri by Tidal Stripping of a Dwarf Galaxy 
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ABSTRACT 

We have investigated whether or not a tidal stripping scenario can reproduce 
the observed surface-brightness profile of uj Centauri using A^-body simulations. 
Assuming that the progenitor of uj Centauri is a dwarf elliptical galaxy, we model 
it with a King model with a core radius being the same as that of uj Centauri. 
A dark matter halo of the dwarf is not taken into account. We consider two 
different models of the Milky Way potential: a singular isothermal sphere and a 
three-component model. The progenitor dwarf is expressed as an A^-body system, 
which orbits in the fixed Galactic potential. The dwarf lost more than 90 per cent 
of its mass during the first few pericenter passages. Thereafter, the mass remains 
practically constant. The final surface- density profile is in good agreement with 
the observational data on uj Centauri, if the pericenter distance of the orbit of 
the progenitor dwarf is around 500 pc. This value is within the error bar of the 
current proper motion data on uj Centauri and Galactic parameters. Although 
our simulation is limited to a King-like progenitor dwarf without a dark matter 
halo, it strongly suggests that the current density profile of uj Centauri is nicely 
reproduced by a tidal stripping scenario, in other words, that uj Centauri can 
plausibly be identified with a stripped dwarf elliptical. 

Subject headings: Galaxy: formation — globular clusters: individual (NGC 5139) 
— methods: N-body simulations 
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1. INTRODUCTION 

uj Centauri (NGC 5139) is, to the best of our knowledge, the most massive globular 
cluster in the Milky Way, with a mass around 5 x lO^Mo (Meylan et al. 1995). Furthermore, 
it shows peculiar chemical and dynamical features, such as the wide spread in metallicity 
distribution of member stars (Freeman & Norris 1981) and the difference in kinematics of 
metal-rich and metal-poor populations (Norris et al. 1997). Some member stars show an 
enhancement of s-process elements produced by AGB stars (e.g.. Smith et al. 2000), which 
means that the potential well must have been deep enough to trap the gas ejected from AGB 
stars (Gnedin et al. 2002). A merger model (Icke & Alcaino 1988) can explain the kinematics 
and metallicity spread, but has difficulty in accounting for the presence of s-process elements. 

A tidal stripping scenario (Zinnecker et al. 1988; Freeman 1993) has been proposed to 
explain all observed features simultaneously. In this scenario, the cluster was born as a 
dwarf elliptical galaxy, which was much more massive than the current u Centauri. As it 
sinks to the Galactic center through dynamical friction, it loses most of its mass through 
the tidal stripping by the gravitational potential of the Milky Way. However, if the central 
density is sufficiently high, the core of the dwarf can still survive as a bound system and 
become a relatively large globular cluster, something like uj Centauri. In this scenario, a 
relatively long star formation history is naturally explained because the potential well of the 
progenitor dwarf elliptical was initially much deeper than that of the present-day uj Centauri. 
In addition, since the progenitor dwarf could be formed through hierarchical merging events, 
signatures of past merging events can also be explained. Furthermore, recently. Martini & 
Ho (2004) reported that some of massive globular clusters in NGC 5128 may also be such 
tidally stripped dwarf ellipticals. 

From the viewpoint of dynamics, whether or not the stripping scenario works depends 
on the following two questions. The first one is whether the dynamical friction can bring 
the progenitor dwarf to the current orbit of uo Centauri. The second one is whether the 
tidal stripping can actually produce the observed surface brightness of uo Centauri. The first 
question has been addressed by recent work (Zhao 2002; Tsuchiya, Dinescu, & Korchagin 
2003; Bekki & Freeman 2003). However, the second question has not been studied so far, 
despite its significance to the validity of the stripping scenario. 

In this paper, we have investigated whether a dwarf galaxy with the orbit similar to 
that of the present-day uj Centauri will evolve to have the spatial structure which agrees 
with that of uj Centauri using A'^-body simulations. In § 2, numerical models and method 
are described. Results are shown in § 3. A discussion is given in § 4. 
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2. MODELS AND METHOD 

We consider the dynamical evolution of a dwarf elliptical in a fixed potential. In this 
paper, we neglect the effects of dynamical friction. This neglect should not affect results as 
far as a progenitor dwarf along the present-day orbit of u Centauri is concerned, since our 
main goal is to see whether or not the tidal field of the Milky Way in the present-day orbit 
of uj Centauri can account for its density structure. In reality, dynamical friction must have 
played a certain role in bringing the progenitor dwarf to its current orbit, but that part of 
evolution is not in the scope of the present paper. 

We use the data compiled by Meylan (1987) for the surface- brightness profile of uj 
Centauri. This profile is well-fitted (except in the outermost region) by a King model with 
a non-dimensional central potential Wq — 5.5 and a core radius Rc — 4.6 pc. We choose 
the total mass of this King model to be identical to the mass suggested from observations, 
5.0 X 10^ Mq. Then, we set up the model of the progenitor dwarf to have the same core 
radius and central velocity dispersion, but with a much more extended halo. We choose a 
King model with Wc = 12 as a model for the dwarf. Since some dwarf elliptical galaxies can 
also be fitted by a King model, we beheve this choice is justified. Tsuchiya, Korchagin, & 
Dinescu (2004) have found that a King-like dwarf elliptical cannot evolve to u; Centauri on 
the basis of numerical simulations of King- and Hernquist-like dwarf galaxies embedded in 
a live Milky Way potential. This point will be discussed in § 4.1. 

The dwarf elliptical is represented by 2^° ~ 1 million particles of equal mass. The initial 
total mass is 1.3 x 10^ Mq. We do not include a dark matter halo of the dwarf, which would 
be more extended than the distribution of visible stars. Since the extended halo would be 
stripped out anyway, we believe the presence of a dark matter halo does not change results 
significantly. In most of our simulations, the final cluster still had more than 50,000 particles, 
with the result that the relaxation time is much longer than the duration of the simulation. 

We use two different models for the gravitational potential of the Milky Way, a singular 
isothermal sphere (run std) and a bulge-disk-halo model (run mwd). The circular velocity 
of the std potential is 220 km sec~^. The mwd potential is almost the same as the model 
used in Johnston, Spergel, & Hernquist (1995). This model consists of three components: 
a Miyamoto-Nagai (1975) disk, a Hernquist (1990) bulge, and a logarithmic halo potential. 
The only difference between the model of Johnston et al. (1995) and ours is the core radius of 
the logarithmic potential. We set the core radius to be 14.0 kpc, while Johnston et al. (1995) 
set it to be 12.0 kpc, so that the circular velocity at the solar circle radius, i.e., 8.0 kpc, is 
220 km sec~^. The potentials and parameters used in simulations are listed in Table 1. The 
circular speed curves of these models are plotted in Figure 1. 
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We use the following kinematical data on to Centauri: the distance from the Sun is 
4.9 kpc, the proper motions arc (/i^ cos 5, /i^) = (—5.08 ± 0.35, —3.57 ± 0.34) mas yr~^, and 
the radial velocity is 232.2 ± 0.7 km scc^^ (Dinescu et al. 1999). For the most probable 
values of the proper motions and the radial velocity in the isothermal model (run std), the 
pericenter and apocenter distances are 1.0 kpc and 6.4 kpc, respectively. We made two other 
runs, in which we expand the initial dwarf by a factor of 1.5 (run rl5) and 2 (run r20), to 
see the effect of the change in the tidal force at the pericenter. For the singular isothermal 
model, the only characteristic scale in the experiment is the core radius of the cluster. A 
linear expansion of the cluster by a factor of 1.5 or 2 can therefore equally well correspond 
to a reduction in the pericenter distance of the unmodified cluster by a factor of 1.5^/^ ~ 2 
(run rl5) or 2^/^ ~ 3 (run r20), which corresponds to ~ 600 pc and ~ 400 pc, respectively. 

We use a hierarchical tree algorithm (Barnes & Hut 1986) on the GRAPE-6 hardware 
(Makino et al. 2003) with an opening angle being 0.5. We employ a Plummer softening of 
0.5 pc, which is roughly one-tenth of the observed core radius of a; Centauri, 4.6 pc. The 
equations of motion are integrated using a leap-frog method with a constant time-step of 
At = 1.7 X 10*^ yr. This time-step is comparable to the time required for a particle with 
the maximum circular velocity around uj Centauri to cover the softening length. The total 
energy was conserved to better than 0.005 per cent in all simulations. 



Figure 2 shows snapshot images from run std, along with the orbit of the cluster center 
for the periods of 0.1 Gyr before and after the time for the snapshot (whenever the orbit is 
available). The cluster center is defined as the position of the particle with the minimum 
potential energy for the iV-body particles. We include all particles to calculate the potential. 
Removing unbound particles did not affect the result, even after more than 90 per cent of 
the particles were stripped. We can see that the stripped stars remain close to the cluster 
orbit, and form numerous ripple structures. 

Figure 3 shows the time evolution of the cluster mass. Here, we define the cluster mass 
simply as a mass within 2rt, where rj is the tidal radius at the apocentric position expressed 
as 



Here Rq is the distance from the Galactic center to the cluster, Mq, is the Galaxy mass within 
-Rg) and vrtc is the cluster mass. We assume that the stars are stripped from the cluster if their 
distances from the cluster center are larger than two apocentric tidal radii (e.g., Baumgardt 
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RESULTS 
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& Makino 2003). For the three-component model, we have determined the cluster mass in 
the same way. The cluster mass exhibits periodic sudden drops, which correspond to the 
pericentric passage. Thus, almost all the mass loss occurred at the pericentric passage. Since 
the cluster is expanded in runs rl5 and r20, the mass loss is larger for these runs than that 
for run std. There is very little difference between the result of run std and that of mwd, 
suggesting that the disk does not contribute significantly to the tidal force. This result is of 
course not surprising, because the disk gravity is small at a pericenter distance of 1 kpc. 

Figure 4 shows the final surface-density profiles for all runs. We assume a constant 
mass-to-light ratio, since the estimated relaxation time of uo Ccntauri is too long for the 
mass segregation to be visible (e.g., Gicrsz & Heggie 2003). Clearly, the result of run std 
does not agree with observations. Replacing the Galactic potential with a more realistic 
one (run mwd) does not improve the situation. On the other hand, changing the pericenter 
distance has a drastic effect, and the agreement of run r20 with observations is very good. 
The result of rl5 falls between std and r20. 

We conclude that a tidal stripping scenario involving a King-like progenitor without a 
dark matter halo can reproduce the observed luminosity profile quite well, provided that the 
pericentric distance of uj Centauri is around 400 pc. 

4. DISCUSSION 

4.1. Compcirison with Previous Work 

Our conclusion that a King-like dwarf can evolve to uj Centauri apparently contradicts 
with the result by Tsuchiya et al. (2004) that final clusters started from King-like dwarf 
ellipticals are still more massive than 10^ Mq. This difference is mainly due to the high 
velocity dispersion of the initial dwarf models in their simulations. The central velocity 
dispersion in their model K5 reaches ~ 55 km sec~^, while we set it to be 22 km sec~^. This 
means that the cluster mass inside the tidal radius in their model K5 is 2.5^ ~ 16 times 
larger than the mass in our model, even if the pericenter distance is assumed to be equal. In 
terms of the fraction of total mass lost, their results are actually quite consistent with ours. 

4.2. Pericenter Distance 

We have shown that the observed surface-brightness profile of uj Centauri is nicely repro- 
duced by the tidal stripping scenario if the pericenter distance of the cluster is roughly 400 pc. 
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Here, we discuss whether such a small pericenter distance is compatible with observations. 

We have adopted the following error estimates for observational data. For the errors 
in the proper motion and the radial velocity, we have used the data by Dinescu et al. 
(1999): (U, V, W) = (-64 ± 11, -254 ± 9, 4 ± 10) km sec-^ Here, (U, V, W) is the relative 
velocity of u; Centauri to the local standard of rest. The positive directions of the U, V, 
and W components are, respectively, outward from the Galactic center, toward the Galactic 
rotation, and toward the Galactic north pole. For the Galactic parameters, Reid et al. (1999) 
estimated the ratio of the circular speed to the solar circle radius, 0o/-Roi on the basis of 
the proper motion study of Sagittarius A* by using VLBA, and concluded that the circular 
speed is 219 ± 20 km sec^^, assuming that the solar circle radius is Rq = 8.0 kpc. For the 
error in the circular speed at the solar circle, we have adopted the value of 220 ±20 km sec~^. 
Since the solar circle radius would have an error comparable to that in the circular speed, 
we have used 8.0 ±1.0 kpc for the solar circle radius. 

We calculated the variations of the pericentric distance due to errors in the orbital 
velocity of u; Centauri, the solar circle radius, and the circular speed at the solar circle. The 
results are summarized in Table 2. The error in the circular speed has a strong effect on the 
pericenter distance, since the orbit of cu Centauri is highly eccentric. A 10 per cent variation 
in the circular speed can lead to 50 per cent change in the pericenter distance. The error 
in the V component also has such large magnification. Thus, our result that the pericenter 
distance must be around 400 pc is compatible with observations. 

We set the total mass of u; Centauri to 5.0 x 10^ M©, following the estimate by Meylan 
et al. (1995). However, since this mass estimate is based on a multi-mass King model, the 
total mass is probably overestimated (see Baumgardt et al. 2003). Indeed, in our best-fit 
King model, the central velocity dispersion reaches 22 km sec~^, while the observed value 
is around 17 km sec^^ (Merritt, Meylan, & Mayor 1997). For the observed value of the 
central velocity dispersion, the total mass is around 3 x 10^ Mq. This change of the cluster 
mass reduces the tidal radius by about 15 per cent (r^ oc m^^). In other words, to yield 
the observed tidal radius with a reduced cluster mass, the periastron distance in the model 
would have to be increased by 30 per cent [from equation (1), Rq oc mZ^^'^\. Therefore, if we 
use this revised estimate of the mass of uj Centauri, the most likely value of the pericenter 
distance is around 500 pc. 
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4.3. Long-Term Orbit Evolution in the Three-Component Model 

In section 3, we have found that the quantitative resuh remains unchanged when we 
change the Galactic potential model from the singular isothermal sphere to the more realistic 
three-component model. However, since we integrated our model only for 0.88 Gyr, we could 
have missed some important contributions of the non-sphericity of the potential. 

We integrated the orbit of a point-mass cluster in the three-component potential for 
~ 1Q13.5 Here, such a long integration time is required since chaotic diffusion time could 
be much longer than a Hubble time. We have found that the minimum pericentric radius 
is ~ 0.9 kpc, although the orbit is chaotic. Consequently, the effect of chaotic diffusion is 
neghgible. 

4.4. Globular Clusters as Probes of Galactic Structure 

We have found that the structure of u) Centauri is nicely explained by the tidal stripping 
scenario, for the orbital parameters within the observational errors. The pericenter distance 
is the most important parameter that determines the structure of the simulated uj Centauri. 
In other words, we can put a fairly tight constraint on the pericenter distance of uj Centauri 
from its internal structure. This constraint can be converted to constraints on the Galactic 
potential, such as the circular speed, if we have high-accuracy data for the proper and radial 
motions of uj Centauri. 

We can probably apply a similar technique to other tidally-limited globular clusters, 
though in this case the modeling is somewhat more complex since we cannot ignore thermal 
relaxation. Clusters in highly eccentric orbits are most useful, since a small change in the 
Galactic parameters results in a large change in the pericenter distance. Thus, the dynamical 
simulation of the evolution of globular clusters in a Galactic tidal field, combined with high- 
accuracy proper motion data which will be available via next-generation astrometry projects, 
will provide us with a unique tool to probe the structure of the Milky Way. 

We thank the anonymous referee for useful comments that helped us to improve the 
clarity of the paper. We are grateful to Prof. S. Hozumi for his critical reading of the 
manuscript. All simulations were run at the University of Tokyo. Some part of data analysis 
was made on workstations at Astronomical Data Analysis Center, National Astronomical 
Observatory of Japan (ADAC/NAOJ). 
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Table 1. The Milky Way models. 



isothermal sphere three component 



GMd 



Vc = 220 km sec-i = 1.0 x lO^^M©, a = 6.5 kpc, b = 0.26 kpc 



Mb = 3.4 X lO^^Mo, c = 0.7 kpc 



$h = In ('^^ + 

ri, = 128 km see-^ d = 11 kpe 



Table 2. The pericenter distance of oj Centauri in the isothermal sphere model. 



observable value pericenter 

velocity (C/, V, W) (-64 ± 11, -254 ± 9, 4 ± 10) km sec"^ 1.0 ± 0.3 kpc 

solar circle radius 7.0 kpc 0.9 kpc 

9.0 kpc 1.0 kpc 

circular speed 240 km sec~^ 0.6 kpc 

200 km sec-^ 1.5 kpc 




Fig. 1. — Circular speed curves of our Galactic models. The dashed, dash-dotted, and dotted 
curves give the contributions of the bulge, disk, and halo, respectively, while the solid curve 
gives the total circular speed in the three-component model. The horizontal dotted hue at 
220 km sec~^ is the rotation curve of the isothermal model. 
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Fig. 2. — Snapshot images from run std. Bound particles are drawn in black, while unbound 
ones are in gray. Simulation time is shown at the top-right corner of each box. The orbit of 
the cluster center for periods of 0.1 Gyr before and after each simulation time is also drawn 
in the solid line. 
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Fig. 3. — Time evolution of the cluster mass for all runs. The sohd, dotted, dash-dotted and 
dashed curves represent the results of runs std, mwd, rl5, and r20, respectively. 
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R (pc) 

Fig. 4. — Final surface-density (brightness) profiles of the clusters for all runs. The curves 
have the same meanings as in figure 3. The open circles with error bars show the surface- 
density profile calculated from the observed surface brightness on the assumption of a con- 
stant mass-to-light ratio. The thin solid line corresponds to the initial model. 



